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ABSTRACT 



IhtravasGular adtrasound (IVUS) is a medical imagmg technique diat 
not only provides three-dimensional information about the blood 
vessel lumen and wall, but also directly depicts atherosclerotic plaque 
stracture and morphology. Automatic processing of large data sets of 
rVTJS data represents an inq}ortant challenge due to ultrasound 
speckle and technology artifacts. A new semi-automatic IVUS 
segmentation model, the fast-marching method, based on grayscale 
statistics of the images, is conipared to active contour segm^tation. 
With fast-marching segmentation, the lumen, intima plus plaque 
stracture, and media contours are cortiputed in parallel. Preliminary 
results of this new IVUS segmentation model agree very well with 
vessel wall contours. Moreover, fast-marching segmentation is less 
sensitive to initializatiGQa with average distance between segmentation 
performed with different initializations < 0.8S % and Haussdorf 
distance < 2.6 %. 
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APPARATUS AND METHOD FOR INTRAVASCULAR ULTRASOUND 



IMAGE SEGMENTATION : A FAST-MAWr RTNQ METHOD 



1. Introduction 

Over the past few years» tiie intravascular idtrasound (IVUS) technology has 
become very useful in research and in some clinical applications, particularly for 
studying adierosclerotic disease. IVUS is a catheter based medical imaging 
technique that produces cross-sectional images of blood vessels. It also provides 
quantitative assessm^t of the vascular wall, information about the nature of 
atherosclerosis lesions as well as the plaque shape and size. IVUS adds 
complementary information to angiography when characterizing atherosclerotic , 
plaque. The tomographic nature of intravascular ultrasound also makes 3D 
reconstruction of the vessel wall possible. Furthermore, quantitative 
measurements of ath^osclerosis disease such as plaque volume, intima-media 
thickness, vascular remodeling and lumen area stenosis can be retrieved from 
IVUS. However, the amoimt of data produced by an intravascular ultrasoimd 
acquisition can be quite large, making the analysis of the data long and 
fastidious. Automatic image segmentation should ease the examination of 
intravascular ultrasound series. On the other hand, because IVUS image quality 
remains poor due to speckle noise, imaging artifects and shadowing of parts of 
the vessel wall by plaque components like calcifications, it is necessary to 
develop specific segmentation methods. 

So far, a number of segmentation techniques have been developed for IVUS 
data analysis. A great portion of this work is based on local properties of in^ige 
pixels: gradient based region growing, gradient active surfaces, pixel intensity 
and gradient active contours, sraph search using first and second derivative 
gradient filtCTs combined to local grayscale properties or Sobel-like edge 
operator, and multiagent image segmentation. Texture information has also been 
investigated as a basis for IVUS segmentation. Finally, optimization of a MAP 
estimator modeling ultrasound speckle and contour geometry was considered. 
Still, in 2001, a clinical expert consensus firora the American College of 
Cardiology reported that no IVUS edge detection method had found widespread 
acceptance by clinicians. 

In the present study, usage of a global image characteristic, i.e. the probability 
density functions of flie IVUS wall components, combined to a &st-marching 
segmentation model, a region-b^ed method handling topological contour 
changes, is proposed to address the IVUS segmentation problem. 
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2. Segmentation Model 
2.1 Fast*Marching Method 

The fast-marching method is a particular case of the level-set approach 
introduced by Osher and Sethian to follow interface propagation. The level-set 
model defines the interface as the zero level of a function ^ of higher dimension. 
The ^ function moves in its normal direction according to a speed function F. 
For a 2D interface, the evolution of ^ is given by equation (1) with initial 
surface ^(x,y,t = 0). 

The level-set model can be applied to image segmentation by inteipreting 
image contours as the propagating interface final position. To achieve this, the 
speed fimction should become close to zero when the propagating front meets 
with the desired contour thus making the interface stop at this position and then 
ending the segmentation process. 

Fast-marching consists of the evolution of an interface propagating imder a 
unidirectional speed function. In this case, the evolving contour must be inside 
tihe region to segment (for a positive speed function or outside for a negative 
one). In the fast-marching foramlation, the evolution of the contour is expressed 
in terms of the arrival time T(x,y) of the contour at point (x,y). The T fonction 
satisfies equation (2), stating that the arrival time difference between two 
adjacent pixels increases as the velocity of the contour decreases. 

|vr|F=i. (2) 

The propagation of the interface is done via the construction of the arrival time 
function. The construction algoridmi selects the intoface point having the 
smallest arrival time and calculates the arrival times of its neighbors, and so on 
until the interface has propagated across the whole image. The interface front is 
considered as stationary when the time gradient is sufficiently high. 

Since multiple contours (lumen, intima and media) must be id^tified on the 
rVUS sCTies, image segmentation is done via a multiple interface extension of 
the fast-marching algorithm. A speed function is then defined for each 
propagating interface. Also, every interface associated to the vessel wall 
contours evolves at a velocity defined in terms of the probability density 
function Pi^i of the corresponcfing anatomical structure. The propagation speed 
of interfece m eL, where L is the set {1, 2, N} of the JV evolving interfeces, 
is giv^ by equation (3). 
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Where is the grayscale of pixel s at position (Oc, 3;^ in image /, PiflJ and 
Pfn(Is) are the measured occiuring probabilities of pixel Is in region / and m. 
Because the occurring probability is more significant for a region tilian for a 
single pixel, the speed fimction is calculated over a certain number Ny of 
neigihbors, which are the 8-connected pixels around (x, y). According to equation 
(3), the m inteifece velocity will always be positive and take higher values when 
it is inside a region having a grayscale distribution close to P^. This velocity 
function was chosen because of its general form (no specific PDF is required) 
and because it provides neighborhood averaging. 

To allow comparison, the segmentation model has also been inq>lemented 
using a gradirat based speed fimction of flie form F(x.y)^\/\\'^VI(x,y)\. Notice 
lhat Ihe multiple interface &st-niarching enables simultaneous segmentation of 
different parts of the vessel wall. 



2*2 Active Contours Model 

To con^are the new IVUS segmentation method based on the fast-marching, 
active contours or snakes, introduced by Kass et ah have also been implemoited 
using PDFs and image gradient The snake genoal energy fimction is given by 
equation (4). For PDFs based snakes, the ext^al energy term is given by 
equation (7). 

^d^mm^ ) = - K - V,., || • 




Where Fis the set of n points v/ making the contour, a, jff and y are weighting 
factors, Is is the intensity of pixel Edjnean is an energy term keeping distance 
between nei^boring points close to the mean distance Ecurv is the energy term 
minimizing the curvature. Pin and Pout are the PDFs inside and outside the 
contour. Eext is the energy term related to image Fs characteristics- The energy 
fimction is calculated over regions vi" and vf^^. These regions are delimited by 
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the current point, the two closest contour points and these three points extended 
by five pixels inside and outside the contour along its normal direction. For the 
gradient implementation of snakes, flie external energy function is given by 
equation (8), 



Where g is the gradient, g,nin and g^ax are tihie minimmn and maximum values of 
gradient among the v/ nei^bors. Segmentation is jSnished, when minimal energy 
is reached. The energy function is nadnimized, by evaluating iteratively the 
energy for three contour point neighbors in normal direction, inside and outside 
the snake. Neighbora minimizing EsnakeCV) are chosen to replace contour points 
of higher energy. Segmentations of the lumen, intima plus plaque structure and 
niedia are performed one at a time. 

2.3 Data Initialization 

The fast-marching segmentation requires an initialization which was done by the 
user on longitudinal images of the IVUS series. Contours were manually traced 
on 3 to 4 different cut planes, creating a set of contour points on each IVXJS 
frame. Then, splines passing tiirough these points were computed and used as 
initial contours, directly for the active contours model and slightly shrunk for the 
fast-marching algoriflun (because the propagating interface must be inside the 
region to segment). 

2A Probability Density Function Estimation 

For intravascular ultrasound images, the limien, intima plus plaque stmcture, 
and media PDFs can be naodelcd by a shifted Rayleigh law. The PDF, with 
parameters {miris a}, is given by equation (9). 



with u > mifiy a > 0^ and flie variance = a^{4-rr)f2. By measuring the standard 
deviation and minimum gray level value of each area outlined by the user's 
initialization, the parameters {min, a} were estimated for each region of interest 

2.S Materials 

The active contours and the multiple interface fast-marching methods were used 
to segment the lumen, intima plus atherosclerotic plaque stracture, and media of 



(8) 




(9) 
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a 200 rVUS image series of a femoral artery. Data were acquired with a Jomed 
equipment (In-vision gold, Helsingborg Sweden) using a 20 MHz transducer. 
Images of size 10 x 10 mm were digitized on 384 x 384 pixel matrices and 
stored using the Dicom standard. The acquisition was done at a 10 images/sec 
frame rate and the catheter pullback velocity was set to 1 mm/sec generating 0.1 
mm thick slices. Image acquisition was not ECG-gated. 

3. Results 

In most cases, a quaUtative analysis of the PDFs based fast-marching 
segmentation revealed detected contours that were very close to the Ixraien, 
intima and media boundaries. A typical segmentation result is shown in Fig. 1 . ' 




Fig. 1. (a) Typical IVUS cross-sectional image and (b) IVUS image with PDFs 
based fast-marching detected contours of lumen, intima plus plaque and media. 

Fig. 2 and Fig. 3 show a volumic representation of the lumen and media 
contours also obtained with PDFs based fast-marching segmentation. On each 
figure, 20 detected contours are drawn over two longitudinal cut planes of the 
200 image series. The other methods (fast-marching with gradient and active 
contour with both gradient and PDFs) also provided similar qualitative results 
(data not shown). " 

To quantify the variability of the segmentation under different initiahzations, 
fast-marching and active contour methods, both based on gradient and PDFs, 
were applied to a 15-image sequence with Ifaree different sets of initial contours. 
Table 1 shows the average and Haussdorf distances (maximum distance to the 
closest point) between resulting contours. The distances are presented in 
percentage of the image size. 




to'i^T^?''*'''^' fast-marching segmentation displayed 

over two IVUS images and two longitudinal cuts. 




t^^^TQ^'^''^^'^' ^^'"^ fast-marching segmentation displayed 

over two rvUS images and two longitudinal cuts. 

I>h n/f^'^T-!;'?'^'^'^ Haussdorf distance (HD) between contours 

with different mitializations in percent. 
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4. Discussion 

Tlie goal of this , work was to demonstrate the IVUS segmentation potential of 
the fast-marching method, and the usefuhiess of region statistics such as 
probabihty density functions of vessel wall anatomical structures. Fig. 1 showed 
tiiat vessel wall boundaries can be identified even when the contrast is not very 
high and when the shape is irregular as for the luminal contour. Figs 2 and 3 
demonstrated that segmentation results followed the desired boundaries for a 
whole volume. In the case of the snake model, the boundaries seemed to be 
more affected when they were not well defined (results not shown). 
QuaHtative analysis of the detected contours indicates that the fast-marching is 
an accurate segmentation method for intravascular ultrasound imaging, but this 
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Qualitative analysis of the detected contours indicates that the fest-marching is 
an accurate segmentation method for intravascular ultrasound imaging, but this 
should be further confirmed from validation with several manually traced 
contours by independent e^qperts. This step would also enable quantitative, 
accuracy comparison of active contour and evolving interface algorithms. 
Nevertheless, quantitative evaluation and conq)arison of the variability or 
precision of active contours and fast-marching under different initializations was 
perfomied in Table L The average and Haussdorf distance comparison metrics 
were chosen instead of area or perimeter diflferraces because fhey directly depict 
point to point contour variations. Table 1 indicates that PDFs based fast- 
marching has the smallest Haussdorf distance, which remained under 2.6% for 
all boundaries, with a relatively small average distance between contours, of less 
than 0.84%. Thus, fast-marching detected contours had small variations when 
initialized differently and the maximum distance to the closest point, 
representing the worst case, stayed low. PDFs active contours also had a small 
average distance of less than 0.78%, but the worst variation could be higher than 
for the evolving interfaces. Table 1 also showed that gradient based 
segmentations, independently of the method, produced more variations (less 
precision) than the ones using grayscale PDFs. This indicates that the PDF 
information seems more constant and possibly more reliable over the whole 
image than the gradient. 

These preliminary segmentation results showed that the fast-marching approach 
is a promising technique for IVUS image processing, in the near future, a need 
would be to automate the initialization step to minimize user's interactions by 
fitting a mixture of PDFs directly to the IVUS image, and then by roughly 
estimating an initial segmentation. This way it is hoped that the PDFs based fast- 
marching algorithm may become a reliable fully automatic segmentation tool of 
choice for IVUS, 
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CLAIM; 

1 . Method for intravascular ultrasoimd image segmentation^ characterized by 
the &ct that it comprises a semi-automatic intravascular ultrasound 
segmentation model based on grayscale statistics of images. 
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